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ABSTRACT 

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high 
sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, 
along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical 
polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile 
of the B21 1 filament approaches a power-law behavior p oc r - 2 0±0 - 4 a t large radii and that the temperature profile exhibits a marked drop at small 
radii. The observed density and temperature profiles of the B21 1 filament are in good agreement with a theoretical model of a cylindrical filament 
undergoing gravitational contraction with a polytropic equation of state: P oc p y and T oc p? _1 , with y=0.97±0.01<l (i.e. not strictly isothermal). 
The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further 
suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in 
this picture are -0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations. 
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1. Introduction 

A growing body of evidence indicates that interstellar filaments 
play a fundamental role in the star formation process. In particu- 
lar, the results from the Herschel Gould Belt survey (HGBS) 
confirm the omnipresence of parsec-scale filaments in nearby 
molecular clouds and suggest that the observed filamentary 
structure is directly related to the formation of prestellar cores 
(Andre et al. 2010). While molecular clouds such as Taurus were 
already known to exhibit large-scale filamentary structures long 
before Herschel (cf. Schneider & Elmegreen 1979; Goldsmith 
et al. 2008), the Herschel observations now demonstrate that 
filaments are truly ubiquitous in the cold interstellar medium 
(ISM) (see Men'shchikov et al. 2010; Molinari et al. 2010; 
Arzoumanian et al. 201 1). Furthermore, the Herschel results in- 
dicate that the inner width of the filaments is quasi-universal 
at ~0.1 pc (Arzoumanian et al. 2011). The characteristic fila- 
ment width corresponds to within a factor of ~ 2 to the sonic 
scale around which the transition between supersonic and sub- 
sonic turbulent motions occurs in diffuse, non-star-forming gas 
and a change in the slope of the linewidth-size relation is ob- 
served (cf. Goodman et al. 1998; Falgarone et al. 2009; Federrath 
et al. 2010). This similarity suggests that the formation of fila- 
ments may result from turbulent compression of interstellar gas 
in low-velocity shocks (cf. Padoan et al. 2001). Alternatively, the 
characteristic width may also be understood if interstellar fila- 
ments are formed as quasi-equilibrium structures in pressure bal- 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



ance with a typical ambient ISM pressure P ext ~2-5xl0 Kcm 
(Fischera & Martin 2012, Inutsuka et al. in prep.). 

The HGBS observations also show that the prestellar cores 
identified with Herschel in active star-forming regions such as 
the Aquila Rift cloud (cf. Konyves et al. 2010) are primarily 
located within the densest filaments for which the mass per 
unit length exceeds the critical value (e.g., Inutsuka & Miyama 
1997), M line , crit =2 c 2 JG~ 15M /pc, where c s ~0.2 km/s is the 
isothermal sound speed for T~1Q K. These Herschel results 
support a scenario according to which core formation occurs 
in two main steps (e.g., Andre et al. 2010, 2011). First, large- 
scale magneto-hydrodynamic (MHD) turbulence gives rise to a 
web-like network of filaments in the ISM. In a second step, grav- 
ity takes over and fragments the densest filaments into prestel- 
lar cores via gravitational instability. Indirect arguments fur- 
ther suggest that dense, self-gravitating filaments, which are ex- 
pected to undergo radial contraction (e.g. Inutsuka & Miyama 
1997), can maintain a constant central width of 0.1 pc if they ac- 
crete additional mass from their surroundings while contracting 
(Arzoumanian et al. 201 1). 

In this letter, we present new Herschel observations taken 
as part of the HGBS toward and around the B211/B213 fila- 
ment in the Taurus molecular cloud, which suggest that this fil- 
ament is indeed gaining mass from a neighboring network of 
lower-density striations elongated parallel to the magnetic field 
(see also Goldsmith et al. 2008). Owing to its close distance to 
the Sun (<f~140 pc - Elias 1978), the Taurus cloud has been 
the subject of numerous observational and theoretical studies. 
In particular, it has long been considered as a prototypical re- 
gion and has inspired magnetically-regulated models of low- 
mass, dispersed star formation (e.g., Shu et al. 1987; Nakamura 
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& Li 2008). As pointed out by Hartmann (2002), most of the 
young stars in Taurus are located in two or three nearly parallel, 
elongated bands, which are themselves closely associated with 
prominent gas filaments (e.g., Schneider & Elmegreen 1979). 
The B211/B213 filament discussed here corresponds to one 
of these well-known star-forming filaments (see also Schmalzl 
et al. 2010; Li & Goldsmith 2012). 



2. Herschel observations and data reduction 

The B211/B213+L1495 area (~6°x2.5°) was observed with 
Herschel (Pilbratt et al. 2010) as part of the HGBS in the Taurus 
molecular cloud (see Kirk et al. 2012, for a presentation of 
early results obtained towards two other Taurus fields). Each 
field was mapped in two orthogonal scan directions at 60" s _1 , 
with both PACS (Poglitsch et al. 2010) at 70 fim & 160 fim and 
SPIRE (Griffin et al. 2010) at 250 /mi, 350 /zm & 500 /mi, us- 
ing the parallel-mode of Herschel. For B211+L1495, the North- 
South scan direction was split into two observations taken on 12 
February 2010 and 7 August 2010, while the East-West cross- 
scan direction was observed in a single run on 8 August 2010. 
An additional PACS observation was taken on 20 March 2012 in 
the orthogonal scan direction at 60" s to fill a gap found in the 
previous PACS data. 

The PACS data reduction was performed in two steps. The 
raw data were first processed up to level-1 in HIPE 8.0.3384, 
using standard steps in the pipeline. These level-1 data were then 
post-processed with Scanamorphos version 16 (Roussel 2012), 
to remove glitches, thermal drifts, uncorrected 1/f noise, and 
produce the final maps. The SPIRE data were reduced with HIPE 
7.0.1956 using the destriper-module with a linear baseline. The 
final SPIRE maps were combined using the "naive" map-making 
method, including the turn-around data. 

Zero-level offsets were added to the Herschel maps based on 
a cross-correlation of the Herschel data with IRAS and Planck 
data at comparable wavelengths (cf. Bernard et al. 2010). The 
offset values were 2.3, 37.3, 40.2, 24.8, and 11.5 MJy/sr at 70, 
160, 250, 350, and 500 /mi, respectively. A dust temperature map 
and a column density map (rd ust and Nn 2 - see online Fig. 1) 
were then derived from the resulting images at the longest four 
Herschel wavelengths. The Nn 2 map was reconstructed at the 
18.2" (0.012 pc at 140 pc) resolution of the SPIRE 250 //m data 
using the method described in Appendix A. 

3. Analysis of the filamentary structure 

In order to facilitate the visualization of individual filaments, 
we performed a "morphological component analysis" (MCA) 
decomposition of the Herschel column density map on a basis 
of curvelets and wavelets (e.g., Starck et al. 2003, 2004). The 
decomposition was made on 6 scales and 100 iterations were 
used 1 . The map corresponding to the sum of all 6 curvelet com- 
ponents, shown in Fig. 2, provides a high-contrast view of the 
filaments after subtraction of the non-filamentary background 
and most of the compact cores. Figure. 2 shows that the B211 
filament is surrounded by a large number of lower-density fila- 
ments or striations which are oriented roughly perpendicular to 
the main filament. It is important to stress that the striations can 
also be seen in the original maps (cf. Figs, la and 3a) and that 



1 The MCA software employed for this decomposition 
is publicly avaialble from the "Inpainting routines" link on 
http://irfu.cea.fr/Phocea/Vie_des_labos/Ast/ast_visu.php?id_ast=1800. 
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Fig. 2. High-resolution (18.2") column density map of the Taurus 
B211+L1495 field derived from the Herschel data. The contrast of fil- 
amentary structures has been enhanced using a curvelet transform (cf. 
Starck et al. 2003). Given the typical width ~0.1pc of the filaments 
(Arzoumanian et al. 2011), this map is approximately equivalent to a 
map of mass per unit length along the filaments. The color bar on the 
right shows a line-mass scale in units of the thermal critical line mass of 
Inutsuka & Miyama (1997), which we estimate to be accurate to better 
than a factor of ~ 2 according to a detailed analysis of the radial profiles 
of the filaments (Arzoumanian et al., in prep.). Note that the main B21 1 
filament is thermally supercritical, while the mass per unit length of the 
faint striations is an order of magnitude below the critical value. 

the curvelet transform was merely used to enhance their con- 
trast. Some of these striations are also visible in the 12 CO(1-0) 
and 13 CO(1-0) maps of Goldsmith et al. (2008) at 45" resolu- 
tion. Following Andre et al. (2010), the column density map of 
Fig. 2 was also converted to an approximate map of mass per 
unit length along the filaments by multiplying the local column 
density by the characteristic filament width of 0. 1 pc (see color 
scale on the right of Fig. 2). It can be seen in Fig. 2 that the mass 
per unit length of the main B2 1 1 filament exceeds the thermal 
value of the critical line mass, Mii ne cr i t =2Cs/G, while the mass 
per unit length of the striations is an order of magnitude below 
the critical value. Assuming that the non-thermal component of 
the velocity dispersion inside the filaments is small compared to 
the sound speed, as suggested by the results of millimeter line 
observations (see Hacar & Tafalla 2011, for the L1517 filament 
and Arzoumanian et al. 2012) we tentatively conclude that the 
B21 1 filament is globally gravitationally unstable, while the per- 
pendicular striations are not. 

3.1. A bimodal distribution of filament orientations 

To analyze the distribution of filament orientations in a quantita- 
tive manner, we applied the DisPerSE algorithm (Sousbie 201 1) 
to the original column density map (Fig. la), in order to pro- 
duce a census of filaments and to trace the locations of their 
crests. DisPerSE is a general method based on principles of com- 
putational topology and it has already been used successfully 
to trace the filamentary structure in Herschel images of star- 



2 



P. Palmeirim et al.: Herschel view of the Taurus B21 1/3 filament and striations: Evidence of filamentary growth? 




Right Ascension (J2000) Right Ascension (J2000) 

Fig. 3. (a) Herschel/SPIKE 250 jim image of the B21 1/B213/L1495 region in Taurus. The light blue and purple curves show the crests of the 
B213 and B21 1 segments of the whole filament discussed in this paper, respectively, (b) Display of optical and infrared polarization vectors from 
Heyer et al. (2008), Heiles (2000), and Chapman et al. (2011) tracing the magnetic field orientation in the B211/L1495 region, overlaid on our 
Herschel/SPIRE 250 yum image. The plane-of-the-sky projection of the magnetic field appears to be oriented perpendicular to the B211/B213 
filament and roughly aligned with the general direction of the striations overlaid in blue. The green, blue, and black segments in the lower right 
corner represent the average position angles of the polarization vectors, low-density striations, and B21 1 filament, respectively. 




50 100 150 

Position Angle 



Fig. 4. Histogram of orientations for the low-density striations identi- 
fied with DisPerSE in the B211+L1495 field (displayed in blue). The 
position-angle distribution of available optical polarization (Heyer et al. 
2008, Heiles 2000) and infrared vectors (Chapman et al. 201 1) are also 
shown (green dashed histogram). Gaussian fits to these distributions are 
superimposed, indicating a peak position angle of 34°±13° for the stri- 
ations and 26°±18° for the B-field polarization vectors. The B211 fil- 
ament has a mean position angle of 11 8° ±20° (black triangle and hor- 
izontal error bar) and is thus roughly perpendicular to both the low- 
density striations and the local direction of the magnetic field. 



forming clouds (e.g., Arzoumanian et al. 2011; Hill et al. 2011; 
Peretto et al. 2012; Schneider et al. 2012). Using DisPerSE with 
a relative 'persistence' threshold of 10 21 cirT 2 (~ 5<x in the map 
- see Sousbie 2011 for the formal definition of 'persistence') 
and an absolute column density threshold of 1-2 xlO 21 cm -2 , we 



could trace the crests of the B21 1 filament and 44 lower density 
filamentary structures (see Fig. 3). Due to differing background 
levels on either side of the B211/3 filament (see Fig. la), we 
adopted different column density thresholds on the north-eastern 
side (2xl0 21 crrT 2 ) and south-western side (10 21 cm -2 ). The re- 
sults of DisPerSE were also visually inspected in both the origi- 
nal and the curvelet column density map, and a few doubtful fea- 
tures discarded. The mean orientation or position angle of each 
filament was then calculated from its crest (see Appendix A of 
Peretto et al. 2012 for details). Figure 4 shows the resulting his- 
togram of position angles. In this histogram, the low-density stri- 
ations are concentrated near a position angle of 34°+ 13°, which 
is almost orthogonal to the B211 filament (PA. =1 18°+20°). 
Interestingly, the position-angle distribution of available opti- 
cal polarization vectors (Heiles 2000; Heyer et al. 2008), which 
trace the local direction of the magnetic field projected onto the 
plane-of-sky, is centered on PA. = 26°±18° and thus very sim- 
ilar to the orientation distribution of the low-density striations 
(see Fig. 4). Figure 3b further illustrates that the low-density stri- 
ations are roughly parallel to the B-field polarization vectors and 
perpendicular to the B21 1 filament. 

3.2. Density and temperature structure of the B21 1 filament 

Using the original column density and temperature maps derived 
from the Herschel data (see Appendix A), we produced radial 
column density and temperature profiles for the B211 filament, 
following the same procedure as Arzoumanian et al. (2011) for 
IC5146. We first determined the direction of the local tangent 
for each pixel along the crest of the B21 1 filament as traced by 
DisPerSE. For each pixel, we then derived one temperature pro- 
file and one column density profile in the direction perpendicu- 
lar to the local tangent. Finally, by averaging all individual cuts 



3 



P. Palmeirim et al.: Herschel view of the Taurus B21 1/3 filament and striations: Evidence of filamentary growth? 



15 




Radius [pc] Radius [pc] 

Fig. 5. (a) Mean radial column density profile observed perpendicular to the B21 1 filament and displayed in log-log format, for both the Northern 
(blue curve) and the Southern part (red curve) of the filament. The yellow area shows the (±1<t) dispersion of the distribution of radial profiles 
along the filament. The inner solid purple curve shows the effective 18.2" HPBW resolution (0.012 pc at 140 pc) of the column density map (online 
Fig. 1 - see Appendix A for details) used to construct the profile. The northern and southern column density profiles are very similar up to r~0.4 pc 
(vertical dashed line) and differ significantly only for r>0.4pc, due to different background levels on either side of the filament. The dashed 
black curve shows the best-fit Plummer model (convolved with the 18.2" beam) described by Eq. (1) with p=2. 0+0.4 and ^^=0.03+0.01 pc for 
r<0.4pc, and including a separate linear baseline on each side representing the background for r>0.4pc [see Eq. (Bl) in Appendix B for details]. 
The dashed curve in light green shows a Gaussian fit to the central part of the profile (mean deconvolved FWHM width ~0.09±0.02 pc). (b) Mean 
dust temperature profile measured perpendicular to the B21 1 filament and displayed using a linear scale (black curve). The solid red curve shows 
the best model temperature profile obtained by assuming that the filament has a density profile given by the Plummer model shown in (a) and 
obeys a polytropic equation of state, P oc p y p , and thus T(r) oc p p (r) (y ~ I) . This best fit has y=0.97±0.01. 



along the crest, we obtained a mean column density and a mean 
temperature profile for the B211 filament (Fig. 5). 

To characterize the resulting column density profile we made 
use of an analytical model of an idealized cylindrical filament. 
This model features a dense, flat inner portion and approaches a 
power-law behaviour at large radii. Analytically, it is described 
by a Plummer-like function of the form (cf. Nutter et al. 2008; 
Arzoumanian et al. 201 1) : 

where p c is the central density of the filament, R^ at is the radius 
of the flat inner region, p is the power-law exponent at large radii 
(r»Rfl at ), A p = x B J is a finite constant factor (for 

p>l) that takes into account the filament's inclination angle to 
the plane of the sky (here assumed to be /=0°), and B represents 
the Euler beta function (cf. Casali 1986). The density structure 
of an isothermal gas cylinder in hydrostatic equilibrium follows 
Eq. (1) with p = 4 (Ostriker 1964). 

According to the best-fit model of B21 1 (cf. Fig. 5a), the di- 
ameter of the flat inner portion is 2/?fl at =0.07±0.02pc, which 
is well resolved compared to the 0.012 pc (or 18.2") resolu- 
tion of the column density map. The power-law regime at large 
radii is p oc r -2.o±o.4^ w hich is significantly shallower than the 
steep p oc r -4 profile expected for unmagnetized isothermal 
filaments but would be consistent with models of isothermal 
equilibrium filaments threaded by helical magnetic fields (Fiege 
& Pudritz 2000). Note that our results for the density profile 
of the B211 filament (e.g. mean deconvolved FWHM width 
~0.09±0.02pc) agree with the characteristic width ~0.1 pc 
found by Arzoumanian et al. (2011) for filaments in IC5146, 
Aquila, Polaris, and very similar to the findings of Malinen 



et al. (2012) for the TMC-1 filament (also known as the Bull's 
tail - Nutter et al. 2008). Malinen et al. derived p^2.3 and 
2/?fl a t~0.09pc, also based on Herschel data from the HGBS. In 
addition, the dust temperature profile of the B21 1 filament shows 
a pronounced temperature drop toward the center (Fig. 5b), 
which suggests that the gas is not strictly isothermal 2 . 

A reasonably good model for the structure of the B211 fil- 
ament is obtained by considering similarity solutions for the 
collapse of an infinite cylinder obeying a polytropic (non- 
isothermal) equation of state of the form P oc pi with y%\. 
Kawachi & Hanawa (1998) have shown that the outer density 
profile of such a collapsing cylinder approaches the power law 

2_ 

p oc r 2 -y . For y values close to unity, the model column density 
profile thus approaches p oc r~ 2 at large radii, which is consis- 
tent with the observed profile of the B211 filament. The best 
Plummer model derived above for the density profile [p p (r) - 
Fig. 5a] can be used to estimate the y value which leads to the 
best fit to the observed dust temperature profile 2 under the poly- 
tropic assumption [T(r) oc p p (r) (r_1) ]. The resulting model fit, 
overlaid in red in Fig. 5b, has y=0.97±0.01<l, corresponding to 
p oc r -i-96±o.02 at ] al -g e ra( jjj f or a self-similar, not strictly isother- 
mal collapsing cylinder. 

4. Discussion: Contraction and accretion in B211? 

The results presented in this paper reveal the density and tem- 
perature structure of the Taurus B211 filament with unprece- 



2 We make the approximation that T sas (r) ss Tj us ,(r), which should 
be correct in the inner part of the B21 1 filament at least, since the gas 
and dust temperatures are expected to be well coupled in high-density 
(>3xl0 4 cirr 3 ) regions (see Galli et al. 2002) and the central density of 
the filament is estimated to be « t .«:4.5xl0 4 cirr 3 . 
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dented detail. The shape of the column density profile derived 
for the B211 filament, with a well-defined power-law regime at 
large radii (see Fig. 5a), and the large column density contrast 
over the surrounding background (a factor ~ 10-20, implying a 
density contrast ~ 100-400) strongly suggest that the main fila- 
ment has undergone gravitational contraction. This is also con- 
sistent with the supercritical mass per unit length measured for 
the B211 filament (Mi; ne a:54M0/pc), which suggests that the 
filament is unstable to both global radial contraction and frag- 
mentation into cores (e.g., Inutsuka & Miyama 1997; Pon et al. 

2011) . Observations confirm that the B211 filament has indeed 
fragmented, leading to the formation of several prestellar cores 
(e.g., Onishi et al. 2002) and protostars (e.g., Motte & Andre 
200 1 ; Rebull et al. 20 1 0) along its length. 

The orientation alignment of the striations with optical po- 
larization vectors suggests that the magnetic field plays an im- 
portant role in shaping the morphology of the filamentary struc- 
ture in this part of Taurus. Earlier studies, using similar polariza- 
tion observations of background stars, already pointed out that 
the structure of the Taurus cloud was strongly correlated with 
the morphology of the ambient magnetic field (e.g. Heyer et al. 
2008; Chapman et al. 2011 and references therein). Using the 
Chandrasekhar-Fermi method, Chapman et al. (201 1) estimated 
a magnetic field strength of ~25^/G in the B211 area and con- 
cluded that the region corresponding to the striations seen here 
in, e.g., Fig. 2 and Fig. 3 was magnetically subcritical. 

Theoretical arguments (e.g., Nagai et al. 1998) predict that, 
in the presence of a "strong" magnetic field, low-density, ther- 
mally subcritical filaments such as the striations observed in 
Taurus should be preferentially oriented parallel to the field 
lines, while high-density, self-gravitating filaments should be 
preferentially oriented perpendicular to the field lines. This dif- 
ference arises because low-density structures which are not held 
by gravity have a tendency to expand and disperse, while self- 
gravitating structures have a tendency to contract. In the pres- 
ence of a magnetic field, motions of slightly ionized gas do not 
encounter any resistance along the field lines but encounter sig- 
nificant resistance perpendicular to the field lines. Consequently, 
an initial perturbation in a low-density part of the cloud will tend 
to expand along the field lines and form an elongated structure or 
a subcritical "filament" parallel to the field. Conversely, a self- 
gravitating structure will tend to contract along the field lines, 
forming a condensed, self-gravitating sheet (cf. Nakamura & Li 
2008) which can itself fragment into several supercritical fila- 
ments oriented perpendicular to the field (e.g. Nagai et al. 1998). 
These simple arguments may explain the distribution of filament 
orientations with two orthogonal groups found in Sect. 3.1 (see 
Fig. 4). Other regions imaged with Herschel where a similar dis- 
tribution of filament orientations is observed and a similar mech- 
anism may be at work include the Pipe nebula (Peretto et al. 

2012) , the Musca cloud (Cox et al. in prep.) and the DR21 ridge 
in Cygnus X (Schneider et al. 2010; Hennemann et al. 2012). 

The morphology of the region with a number of low-density 
striations parallel to the magnetic field lines, some of them ap- 
proaching the B211 filament from the side and apparently con- 
nected to it (see Fig. 2), is also suggestive of mass accretion 
along the field lines into the main filament. To test this hypothe- 
sis, we assume cylindrical geometry and use the observed mass 
per unit length M/,„ e to estimate the gravitational acceleration 
g(R)-2GMu„ e {r<R)IR of a piece of gas in free-fall toward the 
B211 filament, where R represents radius. The free-fall veloc- 
ity Vff of gas initially at rest at a cylindrical radius R m it~2 pc 
(corresponding to the most distant striations) is estimated to 
reach Vff=2 [GMn ne ln(R m i t /R)] 1 1.1 km/s when the material 
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Fig. 6. CO emission observed toward and around the B211 filament 
(Goldsmith et al. 2008). Redshifted 12 CO(1-0) emission integrated 
from V , iSj; =6.6km/s to V LSR =1 .4km/s is displayed in red and mostly 
seen to the north-east of the B211 filament. Blueshifted 12 CO(1-0) 
emission integrated from V L 5 R =4.2km/s to V L sr=5.5 km/s is displayed 
in blue and mostly seen to the south-west of the filament. The main body 
of the B211 filament, displayed in green, corresponds to the I3 CO(l- 
0) emission detected between V L 5 fi =5.6km/s to V i- s K =6.4km/s. Both 
12 CO(1-0) and 13 CO(1-0) maps have a spectral resolution of ~0.2 km/s. 

reaches the outer radius R~QA pc of the B211 filament. This 
free-fall estimate is an upper limit since it neglects any form of 
support against gravity. A more conservative estimate can be ob- 
tained by considering the similarity solution found by Kawachi 
& Hanawa (1998) for the gravitational collapse of a cylindri- 
cal filament supported by a polytropic pressure gradient with 
y&l. In this model, the radial infall velocity in the outer parts of 
the collapsing filament is expected to be Vm ~ 0.6-1 km/s when 
y =0.9-0.999 and the gas temperature is ~10K (see Figs. 4 and 
6 of Kawachi & Hanawa). The above two velocity estimates 
can be compared with the kinematical constraints provided by 
the 12 CO(1-0) observations of Goldsmith et al. (2008). It can 
be seen in online Fig. 6 that there is an average velocity differ- 
ence of ~ 1 km/s between the red-shifted CO emission observed 
at Vlsr~1 km/s to the north-east and the B211 filament which 
has Visr~6 km/s. Likewise, there is an average difference of ~1 
km/s between the blue-shifted CO emission observed at Vls«~5 
km/s to the south-west and the B211 filament. Although projec- 
tion effects may somewhat increase the magnitude of the intrin- 
sic velocity difference, we conclude that there is good qualita- 
tive agreement between the estimated inflow velocity in the stri- 
ations and the 12 CO observational constraints. Considering these 
velocities, the current mass accretion rate onto the 4-pc-long fil- 
ament (total mass of ~22OM ) is estimated to be on the order of 
Mi ine = p p (R) x v inf x 27rR w 27-50 M /pc/Myr, where p p (R) cor- 
responds to the density of the best-fit Plummer model at the fila- 
ment outer radius R=0.4pc. This would mean that it would take 
~l-2 Myr for the central filament to form at the current accre- 
tion rate and ~0.8-1.5 Myr for the total mass of the striations 
(~150M o ) to be accreted. The available observational evidence 
therefore lends some credence to the view that the B21 1 filament 
is radially contracting toward its long axis, while at the same 
time accreting additional ambient material through the striations. 
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Appendix A: Derivation of a high-resolution (18.2") 
column density map 

The procedure employed here to construct a column density 
map at the 18.2"resolution of the SPIRE 250 fim data for the 
B211+L1495 region is consistent with, but represents an im- 
provement over, the method used in earlier HGBS papers to de- 
rive column density maps at the 36.3" resolution of SPIRE 500 
fim observations. Following the spirit of a multi-scale decompo- 
sition of the data (cf. Starck et al. 2004), the gas surface den- 
sity distribution of the region, smoothed to the resolution of the 
SPIRE 250 fim observations, may be expressed as a sum of three 
terms: 

S250 = ^500 + (£350 - £500) + (X250 - £350) ■ 

In the above equation, X500, £350, and £250 represent smoothed 
versions of the intrinsic gas surface density distribution Z after 
convolution with the SPIRE beam at 500 fim, 350 fim, and 250 
fim, respectively, i.e.: S 500 = £ * £500, £350 = X * B350, and 

2-250 — ^ * -B250- 

The first term of Eq. (A.l) is simply the surface density 
distribution smoothed to the resolution of the SPIRE 500 fan 
data. An estimate, £500, of this term can be derived from the 
Herschel data using the same procedure as in earlier HGBS pa- 
pers (e.g. Konyves et al. 2010). Briefly, the Herschel images in- 
cluding the zero-level offsets estimated from IRAS and Planck 
(cf. Bernard et al. 2010) are first smoothed to the 500 fim reso- 
lution (36.3") and reprojected onto the same grid. An optically 
thin greybody function of the form I v = B v (T,j) k v £, where I v is 
the observed surface brightness at frequency v, and k v is the dust 
opacity per unit (dust+gas) mass, is then fitted to the spectral 
energy distributions (SEDs) observed with Herschel between 
160 fim and 500 fim, on a pixel-by-pixel basis (four SED data 
points per pixel). This makes it possible to estimate the best- 
fit value £ 5 oo(x,y) and T dy500 (x,y) at each pixel position (x,y). 
The following dust opacity law, very similar to that advocated 
by Hildebrand (1983) at submillimeter wavelengths, is assumed: 
k v = 0.1 X (v/1000 GHzf = 0.1 X (300 fim/Xf cm 2 /g, with 
= 2. 

The second term of Eq. (A.l) may be written as S350 - Z350 * 
G5oo_350, where Gsoo_350 is a circular Gaussian with full width 
at half maximum (FWHM) V36.3 2 - 24.9 2 * 26.4". (To first 
order, the SPIRE beam at 500 fim is a smoothed version of the 
SPIRE beam at 350 fim, i.e., B 500 = B350 * G 50 o_35o-) The sec- 
ond term of Eq. (A.l) may thus be viewed as a term adding in- 
formation on spatial scales accessible to SPIRE observations at 
350/im, but not to SPIRE observations at 500fim. In practice, one 
can derive and estimate £350 of S350 in a manner similar to £500, 
through pixel-by-pixel SED fitting to three Herschel data points 
between 160 fim and 350 fim (i.e., ignoring the lower resolution 
500 fim data point). An estimate of the second term of Eq. (A.l) 
can then be obtained by subtracting a smoothed version of £350 
(i.e., £ 35 o * G 50 o_35o) to £350 itself, i.e., by removing low spatial 
frequency information from 2350. 

Likewise, the third term of Eq. (A.l) may be written as 
£250 - £250 * G 35 o_250, where G 3 50_250 is a circular Gaussian with 
FWHM V24.9 2 - 18.2 2 ~ 17.0", and may be understood as 
a term adding information on spatial scales only accessible to 
Herschel observations at wavelengths < 250/im. In order to de- 
rive an estimate £250 of £250 on the right-hand side of Eq. (A.l), 
we first smooth the PACS 160 fim map to the 18.2" resolu- 
tion of the SPIRE 250 fim map and then derive a color tem- 
perature map between 160 fim and 250/im from the observed 
h5Qfim(x,y)/I\6Qfim(x,y) intensity ratio at each pixel (x,y). The 



SPIRE 250 fim map is converted into a gas surface density map 
(£250), assuming optically thin dust emission at the temperature 
given by the color temperature map and a dust opacity at 250 fim 
K250//m = 0.1 X (300/250) 2 cm 2 /g. An estimate of the third term 
of Eq. (A. 1) can then be obtained by subtracting a smoothed ver- 
sion of £250 (i-e., £250 * G 3 5o_25o) to £ 2 5o_itself, i.e., by removing 
low spatial frequency information from £250- 

Our final estimate £250 of the gas surface density distribution 
at 18.2" resolution is produced by adding up the above estimates 
of the three terms on the right-hand side of Eq. (A.l): 

£250 = £500+(£350 - £350 * G500_35o)+(£250 ~ £250 * G3 5( )_25o) ■ (A 

The resulting 18. 2" -resolution column density map /Vh 2 for 
the B21 1+L1495 region is displayed in online Fig. la in units of 
mean molecules per cm 2 , where £250 = fi m H Nh 2 and H = 2-33 is 
the mean molecular weight. Although this high-resolution map 
is somewhat noisier than its 36.3"-resolution counterpart (cor- 
responding to £500) and has lower signal-to-noise ratio than the 
SPIRE 250 fim image (cf. Fig. 3), due to additional noise com- 
ing from the second and third terms of Eq. (A. 1), the quality and 
dynamic range of the Herschel data are such that the result pro- 
vides a very useful estimate of the column density distribution in 
B211+L1495 with a factor of 2 better resolution than standard 
column density maps derived so far from Herschel observations. 

Because the higher-resolution terms in Eq. (A.l) are derived 
using fewer and fewer SED data points to estimate the effective 
dust temperature for each line of sight, the high-resolution 
column density map is also somewhat less reliable than the 
standard 36.3"-resolution column density map (corresponding 
to £500). To evaluate the reliability of the £350 and £250 maps 
entering the calculation of £250 [cf. Eq. (A. 2)] and derived us- 
ing only three and two SED data points per position, respec- 
tively, we made the following tests in the case of the Taurus 
B211/B213+L1495 data. First, from the Herschel 160 fim to 
350 fim images smoothed to the 500 fim resolution, we derived 
a dust temperature map 7^350 anc * a S a s surface density map 
£^[J using the three Herschel data points between 160 fim and 
350 fim and compared these maps to the maps T^oo and £500 
derived at the same resolution using four SED data points per 
position. In the case of the Taurus data, the T^ 5Q map agrees 
with the rd,5oo map to better than 0.15 K on average (and better 
than 0.8 K everywhere) and £3^ agrees with £500 to better than 
4% on average (and better than 25% everywhere). Likewise, us- 
ing the Herschel 160 fim and 250 fim images smoothed to the 
500 fim resolution, we derived a color temperature map T™ 50 

and a gas surface density map I.™ m the same way as we cal- 
culated £250 above and then compared these maps to the T A 50 o 
and £500 maps. The 7^250 ma P agrees with the 7^500 map to bet- 
ter than 0.15 K on average (and better than 1.2 K everywhere) 
and £^0 agrees with £500 to better than 3% on average (and bet- 
ter than 30% everywhere) 3 . Finally, to test the robustness of the 
18.2"-resolution column density map £250, we smoothed it to 
the 36.3" resolution of the standard column density map and in- 
spected the ratio map between the two, which has a mean value 



3 In pathological situations, such as when warm foreground dust 
emission from a photon-dominated region is present in front of colder 
structures, the difference maps T^f 5Q - T Ai500 and T% f 50 - r di50 o could po- 
tentially be used to improve the estimates of the second and third terms 
of Eq. (A.l). Since these difference maps remain small in the case of 
Taurus, we refrained from using them here and did not apply any cor- 
rection to S350 and S 2 5o- 
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of 1.00 and a standard deviation of 0.04. Within the region cov- 
ered by both PACS and SPIRE, the smoothed version of £250 
agrees with £500 to better than 10%. 

Appendix B: Details of the procedure used to fit the 
column density profile of the B211/B213 filament 

The fitting analysis of the observed column density profile 
(Sect. 3.2 and Fig. 5) was performed using the non-linear least- 
squares fitting IDL procedure MPFIT (Markwardt, C. B. 2008 - 
http://purl.com/net/mpfit). In addition to the Plummer-like cylin- 
drical model filament corresponding to Eq. (1) with three free 
parameters, the background gas was represented by two separate 
linear baselines on either side of the filament. Each side was fit- 
ted independently with the following model function (convolved 
with the beam): 

N° 

X p (r)/nm H = ^ — + Bkg[l]r + Bkg[0], (B.l) 

[l + W/W 2 ] 2 

where N^(- A p p c /? flat / p.m H ), Rfiat, P, Bkg[l], and Bkg[0] were 
treated as five free parameters. Bkg[l] and Bkg[0] are two pa- 
rameters which describe the local background cloud gas. The 
results of these model fits are given in Table B.l for both the 
north-eastern and the south-western side of the B211/B213 fila- 
ment considered as a whole (see Fig. 5a), as well as for the B21 1 
and the B213 segment of the whole filament (see online Fig. 7). 
Due to differing background levels on either side of the filament, 
fitting the two sides of the filament separately gives better results 
than a global fit to both sides simultaneously. 
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Table B.l. Parameters of the Plummer-like model fits to the column density profiles of the B21 1/B213 filament and segments 



Mean 
radial profiles 


Art) (a) 

[10 21 cm" 2 ] 


D 

"flat 

[pc] 


P 


BkgLlJ 

[10 21 cm" 3 ] 


BKgLOJ 
[10 21 cm" 2 ] 


,-2 

X 


north-eastern side of B21 1/B213 
south-western side of B21 1/B213 


14.8+1.9 
14.4+1.4 


0.031 + 0.012 
0.032 + 0.014 


2.03 ± 0.34 
2.00 ± 0.09 


0.86 ± 0.45 
0.0 


0.52 ± 0.80 
0.67 + 0.17 


1.87 
2.12 


north-eastern side of B211 segment 
south-western side of B21 1 segment 


14.9+1.2 
14.0+1.0 


0.050+0.014 
0.057 ± 0.015 


2.51 + 0.45 
2.64 ± 0.50 


0.73 + 0.38 
0.23 + 0.36 


0.95 ± 0.52 
1.31+0.48 


4.37 
3.89 


north-eastern side of B213 segment 
south-western side of B213 segment 


13.0+1.8 
12.7 + 2.0 


0.025 ± 0.008 
0.020 ± 0.007 


1.75 + 0.07 
1.64 + 0.05 


0.95 + 0.11 
0.03 + 0.12 


0.0 
0.0 


1.41 
6.00 



Notes: (a) Mean central column density measured along the crest of the filament (after background subtraction), (b) On the south-western side 
of the B21 1/B213 filament, the background is well described by a constant value, i.e., Bkg[l] = 0. The fitting analysis was performed on the 
mean column density profiles measured on the north-eastern and south-western sides of the global (B211/B213) filament and the two (B211 
and B213) segments (see online Fig. 7). 
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Fig. 1. (a) High-resolution (18.2") column density map of the Taurus B211/B213 region (in units of Nh 2 cnr 2 ) derived from Herschel data as 
explained in Appendix A. (b) Dust temperature map of the Taurus B21 1/B213 region (in K) derived at 36.3" from Herschel data. Comparison of 
the two panels shows how dust temperature and column density are anti-correlated. 
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Fig. 7. Mean radial column density profiles of the B213 (left) and B211 (right) segments of the filament for both the north-eastern (in blue) and 
south-western (in red) sides. In both panels, the black dashed curve represents the mean of the best-fit Plummer models to the north-eastern and 
south-western profiles, truncated at r = 0.4pc where the backgrounds start to diverge significantly between the two sides. The crests defining the 
B213 and B21 1 segments considered here can be seen in Fig. 3a. The B213 segment has a slightly shallower profile (p « 1.7) and a smaller flat 
inner radius (i? flat « 0.025 pc) than the B21 1 segment (p « 2.6 and « 0.06 pc). The detailed parameters of the model fits are given in Table B.l. 



